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ABSTRACT 

The origin of blue straggler stars (BSS) in globular clusters (GCs) is still not fully 
understood: they can form from stellar collisions, or through mass transfer in isolated, 
primordial binaries (PBs). In this paper we use the radial distribution of BSS observed 
in four GCs (M3, 47 Tuc, NGC 6752 and lo Cen) to investigate which formation process 
prevails. We find that both channels co-exist in all the considered GCs. The fraction 
of mass-transfer (collisional) BSS with respect to the total number of BSS is around 
- 0.4 - 0.5 (- 0.5 - 0.6) in M3, 47 Tuc, and NGC 6752. The case of w Cen is 
peculiar with an underproduction of collisional BSS. The relative lack of collisional 
BSS in uj Cen can be understood if mass segregation has not yet driven to the core 
a sizeable number of PBs, which dominate stellar collisions through three- and four- 
body processes. The spatial distribution of BSS provides strong hints to their origin: 
the BSS in the cluster outskirts form almost exclusively from mass transfer in PBs, 
whereas the BSS found close to the cluster core most likely have a collisional origin. 

Key words: stellar dynamics - binaries : general - blue stragglers - globular clusters: 
individual: M3 - globular clusters: individual: 47 Tuc - globular clusters: individual: 
NGC 6752 - globular clusters: individual: lo Cen 



1 INTRODUCTION 

Blue straggler stars (BSS) are stars lying above and blue- 
ward of the turn-off in the color-magnitude diagram (CMD) 
of a star cluster. At least two different processes have been 
proposed to explain their formation (Fusi Pecci et al. 1992; 
Bailyn 1995; Bailyn and Pinsonneault 1995; Procter Sills, 
Bailyn & Demarque 1995; Sills & Bailyn 1999; Sills et 
al. 2000; Hurley et al. 2001). The first, dubbed the mass- 
transfer scenario, suggests that BSS are generated by pri- 
mordial binaries (hereafter PBs) that evolve mainly in iso- 
lation until they start mass-transfer and possibly coalesce 
(McCrea 1964; Carney et al. 2001). The second, known as 
the collisional scenario, states that BSS are the product of 
a merger between two main sequence stars (MSs) in a dy- 
namical interaction that involves a MS-MS collision, most 
likely in a binary-MS encounter (Davies, Benz & Hills 1994; 
Lombard! et al. 2002). The collisional BSS (COL-BSS) differ 
kinematically from the mass-transfer BSS (MT-BSS), since 
they are believed to acquire kicks due to dynamical recoil. 
In both these hypotheses, the resulting BSS have mass ex- 



ceeding the turn-off mass of the cluster and are fueled by 
hydrogen thanks to the mixing of the hydrogen-rich surface 
layers of the two progenitor stars. 

The two aforementioned scenarios do not necessarily 
exclude each other and might co-exist within the same star 
cluster (Ferraro et al. 1993; Davies, Piotto & De Angeh 2004; 
Mapelli et al. 2004) . However, until recently, it has been diffi- 
cult to estimate the relative importance of the two formation 
channels in a given globular cluster (GC). A new approach 
to solve this problem has been presented by Mapelli et al. 
(2004, hereafter paper I), who revised and upgraded an origi- 
nal attempt by Sigurdsson, Davies & Bolte (1994). The basic 
input for the new method is the observation of the shape of 
the BSS radial distribution in a star cluster. After the pio- 
neering works on M3 (Ferraro et al. 1993, 1997a) and M55 
(Zaggia, Piotto & Capaccioli 1997), accurate determinations 
of the spatial distribution of BSS are now being obtained for 
an increasing sample of GCs (Ferraro et al. 2004 for 47 Tuc; 
Sabbi et al. 2004 for NGC 6752; Ferraro et al. 2006a for 
Lu Cen; Warren, Sandquist & Bofie 2006 for M5). A very 
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Figure 1. Comparison between tlie observed surface star density 
profile of the considered GCs (solid line and full circles) and the 
adopted multi-mass King model (dashed line). The units on the 
y-axis are number of stars per arcsec^. The observed profiles are 
taken from Ferraro et al. (1997a) for M3, Ferraro et al. (2006a) 
for uj Cen, Mapelli et al. (2004) for 47 Tuc, Ferraro et al. (2003b) 
for NGC 6752. 



interesting result is emerging from these observations: BSS 
display a clear tendency to follow a bimodal spatial distri- 
bution, with a peak in the core, decreasing at intermediate 
radii and rising again at larger radii. 

In paper I, Mapelli et al. focused on 47 Tuc and showed 
that the bimodal distribution can be reproduced assuming 
a suitable combination of the two proposed mechanisms for 
forming BSS. In particular, it was suggested (paper I) that 
the BSS in the core of 47 Tuc are mainly COL-BSS, while the 
external BSS are MT-BSS. This result agrees with studies 
of the BSS luminosity function in 47 Tuc (Bailyn & Pinson- 
neault 1995; Sills & Bailyn 1999; Sills et al. 2000; Ferraro et 
al. 2003a; Monkman et al. 2006). 

In this paper, we explore two other GCs, whose BSS 
have a bimodal radial distribution. The aim is to test 
whether the conclusions drawn for 47 Tuc can be extended 
to GCs having different mass, concentration, central density, 
and velocity dispersion with respect to 47 Tuc. In particu- 
lar, we study whether the position of a BSS in a GC can 
reliably be regarded as a signature for its origin. We also 
consider the case of lu Cen, whose radial distribution of BSS 
is flat. We follow the same procedure of paper I, carrying out 
dynamical simulations of BSS evolved in the gravitational 
potential well of their associated cluster. Details about the 
simulations are illustrated in Section 2; the sample of the 
investigated clusters is presented in Section 3. The results 
for M3, NGC 6752 and 47 Tuc are discussed in Section 4, 
whilst Section 5 is devoted to the outlier, ui Cen. In Section 
6 we discuss our flndings in comparison with the findings 
of Davies et al. (2004). Finally in Section 7 we present our 
conclusions. 



2 THE SIMULATIONS 

The simulations have been performed with an upgraded ver- 
sion (fully described in paper I) of the code originally de- 
veloped by Sigurdsson & Phinney (1995). This code fol- 
lows the dynamical evolution of a BSS in a static clus- 
ter background, which is modelled using a multi-mass King 
density profile. In this case, once the classes of mass have 
been selected, the cluster background is uniquely deter- 
mined by imposing a central velocity dispersion a, a core 
stellar density ric and a dimensionless central potential Wo 
{Wo = *(0)/{cr)^, where (a) is the mean core velocity dis- 
persion and ^/(O) = 4'(rt) — $(0), with $(r) the gravitational 
potential at the radius r and rt the tidal radius). In practice, 
we have chosen 10 classes of mass (the same classes reported 
in table 1 of paper I) and have adopted the best available 
estimates of a and ric for each of the considered clusters. 
The values of Wo for each GC have been derived by fitting 
the simulated star density profile to the observed one (see 
Section 3 and Table 1 for details). 

Once the background has been determined, the dynam- 
ical evolution of the current population of BSS is simulated 
assuming a value for the ratio between the number of MT- 
BSS and that of COL-BSS. As far as the birth places are con- 
cerned, we assume that COL-BSS are generated exclusively 
in the innermost region, within the core radius (rc), where 
the star density is highest, leading to a high collision rate 
(Leonard 1989; Pooley et al. 2003). MT-BSS can be gener- 
ated everywhere in the cluster, but are expected to be more 
frequent in the peripheral regions, where PBs can evolve in 
isolation, without suffering exchange or ionization by gravi- 
tational encounters (Sigurdsson & Phinney 1993; Ivanova et 
al. 2005). For this reason, in our simulations the MT-BSS, 
formed in PBs, are generated outside the cluster core with 
initial locations distributed in several radial intervals, be- 
tween 1 and 80 rc. Within this region, all initial positions 
are randomly chosen following a fiat probability distribu- 
tion, according to the fact that the number of stars A'^ in a 
King model scales as dN = n(r) dV oc nr'^dr oc dr. 

BSS velocities are randomly generated following the dis- 
tribution illustrated in section 3 of Sigurdsson & Phinney 
(1995, equation 3.3). In addition, we assign a natal kick ve- 
locity Wkick equal to 1 cr to those BSS formed collisionally in 
the core. We tried also different values of Wkick- In agreement 
with paper I, we find that Ukick higher than 2 — 3 o" causes 
the ejection of most of the BSS from the cluster, while there 
are no significant differences in the results for kick velocities 
ranging between and 2 a. The masses of the BSS range be- 
tween 1.2 and 1.5 Mq. In paper I, we extended our analysis 
up to a mass of 2 Mq (Ferraro et al. 1997a; Gilliland et al. 
1998); but now more stringent constraints are coming from 
observations (e.g. the upper limit for the mass of BSS in ui 
Cen is 1.4 Mq; Ferraro et al. 2006a). On the other hand, 
in paper I we found that there is not significant difference 
between the behaviour of a 1.2 and a 2 Mq BSS. Each single 
BSS is evolved for a time ti — fi tiast, where fi is a random 
number uniformly generated in [0, 1] and tiast is the maxi- 
mum lifetime attributed to a BSS. Given the uncertainty on 
the value of tiast, we have performed various sets of runs with 
iiast spanning the range between 1 and 5 Gyr. Consistent 
with the results of paper I, the best fits are all obtained for 
tiast= 1.5 Gyr. For 1 Gyr ^ tiast 2 Gyr the distribution 
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Table 1. Globulcir cluster parameters 





rc (pc) 


(7 (km s ^) 


Uc (stars pc ^) 


Wo 


c 


references** 


NGC 104 (47 Tuc) 


0.47 


10 


2.5 X 10^ 


12 


1.95 


1, 2 


NGC 5272 (M3) 


1.5 


4.8 


6 X 103 


10 


1.77 


1, 3 


NGC 6752 


0.1, O.SS'' 


4.9, 12.4^^ 


2 X 10^ 


13, 12*' 


2.03, 1.95*' 


1, 4 


NGC r,l39 U- Ceil) 


1.1 


17 


r,.o X 10''> 


6.r, 


1.10 





^References: 1 - Dubath et al. (1997); 2- Mapelli et al. (2004); 3- Sigurdsson et al. (1994); 4- Drukier et al. (2003); 

5- Merritt, Meylan & Mayor (1997). 
The double value of Wq and c for NGC 6752 is due to the fact that the profile of this cluster is fit with a double King model. Wq = 13 
and c = 2.03 refer to the inner King (with core radius 5.7"= 0.1 pc for a distance of 4.3 kpc from the Sun; Fcrraro ct al. 2003b); while 

Wo = 12 and c = 1.95 refer to the outer King (with core radius 28"= 0.58 pc; Ferraro et al. 2003b). 
The double value of a for NGC 6752 refers to the two measurements cr=4.9l ^ 4 km s'^ (Dubath et al. 1997) and ct=12.4±0.5 km s"! 

(Drukier et al. 2003). 



does not change dramatically. If we choose a longer lifetime 
(3 Gyr or more), the dynamical friction washes out any peak 
of peripheral BSS; whereas for shorter lifetimes (less than 1 
Gyr) any BSS could be hardly observable today. 

Our code follows the dynamical evolution of the BSS 
in the cluster potential, subject to the action of dynamical 
friction and to the effects of distant encounters (equation 3.4 
of Sigurdsson and Phinney 1995). 

The final positions of the BSS in each run determine 
the simulated radial distribution, to be compared with the 
observed one. We ran 10000 experiments for each of the 
considered cases. Since we studied ~ 10 different cases for 
each cluster, we made about ~ 400000 runs in total. 



3 THE SAMPLE OF GLOBULAR CLUSTERS 

We have studied the distribution of BSS in four different 
GCs: M3, 47 Tuc, NGC 6752 and co Con. They represent a 
very inhomogcncous sample of GCs. M3 is an intermediate 
density (6x 10^ stars pc~^; Sigurdsson et al. 1994), low veloc- 
ity dispersion (4.8 km s~^; Dubath, Meylan & Mayor 1997) 
cluster. 47 Tuc and NGC 6752 are two very concentrated, 
maybe core collapsed clusters, uj Con is the most massive 
among the Galactic GCs (5 x 10*^ M©; Pryor & Meylan 1993; 
Meylan et al. 1995) and it is still a matter of debate whether 
it is a GC or the remnant of a dwarf galaxy (Zinnecker et al. 
1988; Freeman 1993; Ideta & Makino 2004): it is a very loose 
cluster, with a wide metallicity spread and a clear evidence 
for rotation (van Leeuwen et al. 2000; Reijns et al. 2006; van 
de Ven et al. 2006). 

The first step of our work consisted in finding the best 
match between the simulated star density profile of each 
cluster and the observed one. The results arc shown in Fig. I. 
The data of the star density profile of M3, NGC 6752 and 
LO Cen are from Ferraro et al. (1997a, 2003b) and Ferraro et 
al. (2006a), respectively. The star density profile of 47 Tuc is 
the same as published in paper I. Note that we fit the profile 
of NGC 6752 with a double King model, as in Ferraro et al. 
(2003b). The adopted values of a and ric, and the fit values 
of Wo are given in Table I, as well as the concentration 
parameter c (where c = log('rt/rc)) and the fiducial values 
of Tc in parsec. The latter quantities are derived by Eissuming 
a distance from the Sun of 4.6 kpc (47 Tuc, Ferraro et al. 
1999, see also table 3 by Beccari et al. 2006), 10. 1 kpc (M3, 



Fcrraro ot al. 1999), 4.3 kpc (NGC 6752, Ferraro et al. 1999) 
and 5.5 kpc {ui Cen, Bellazzini et al. 2004). 

The data for the radial distribution of BSS are taken 
from Fcrraro et al. (1997a) for M3, Ferraro et al. (2004) 
for 47 Tuc, Sabbi et al. (2004) for NGC 6752 and Ferraro 
et al. (2006a) for u> Cen. As usual, the number of BSS is 
normalized to the number of red giant branch (RGB) stars 
or to the number of horizontal branch (ffB) stars. Even if 
the considered GCs are very different from each other, the 
observed radial distributions of BSS are similar (see Fig. 2), 
with the exception of w Cen. In fact, the radial distribution 
of BSS displays a maximum in the central region of the GC, 
decreases down to a minimum at intermediate radii (of order 
5 — 10 Tc) and increases again at outer radii, ui Cen differs 
as it shows an almost flat distribution. 



4 RESULTS AND DISCUSSION FOR M3, 
47 TUC AND NGC 6752 

We have explored whether the radial distribution of BSS 
observed in our sample of GCs can be reproduced with a 
model which predicts two different classes of progenitors for 
the BSS: (i) PBs evolved in isolation and (it) binary (or 
single) stars which underwent collisions in the high density 
environment of the GC core. Fig. 2 shows the results ob- 
tained for M3, 47 Tuc, NGC 6752 and co Cen. In all cases, 
our adopted model provides a statistically acceptable fit to 
the observations. The best fit value of the fraction ?7mt 
of MT-BSS with respect to the total number of BSS [i.e. 

rjMT = iVMT-BSs/(A^MT-BSS+.^COL-BSs), where iVMT-BSS 

and A^coL-Bss arc the number of mass-transfer and coUi- 
sional BSS, respectively] is reported in Table 2, as well as 
the reduced of the best fit solution. 

In order to determine the best fit solution, we have cre- 
ated a one-dimensional grid of simulations for each cluster, 
allowing the only free parameter in our model^, 77MT, to vary 
between and 1 with an initial step of 0.2. Then, the grid 

^ Other possible parameters (such as the initial velocity of BSS) 
were found to have only minor impact on the results (see Sec- 
tion 2) and have been set to a fixed value in all the considered 
runs. Also, the cluster properties are not free parameters in the 
fitting procedure, because they have been chosen a priori (see Sec- 
tion 3) and kept fbced in all simulations of a given cluster. An ex- 
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Figure 2. Radial distribution of BSS normalized to tlie distribution of horizontal branch (HB) stars (47 Tuc, NGC 6752) or red giant 
branch (RGB) stars (M3, u Cen). Filled circles and solid lines indicate the data taken from Ferraro et al. 1997a (M3), Ferraro ct al. 
2006a (u; Cen), Ferraro et al. 2004 (47 Tuc), Sabbi et al. 2004 (NGC 6752). Open circles and dashed lines indicate the best fits obtained 
from our simulations. For NGC 6752 we adopted rc = 28" (see caption of Table 1). The dotted vertical line in all the panels shows the 
characteristic radius of the zone of avoidance rav (see text and Table 2). In the panel of uj Cen, the open squares and the dot-dashed line 
connecting them report the ratio A^xMM /^RGB , that is the number of X-ray sources detected in lo Cen by XMM (Gendre et al. 2003) 
normalized to the number of RGB stars. 



was iteratively refined down to a step of 0.01 for the values 
surrounding the best fit. 

The values of the reduced were derived by com- 
paring each simulation with the observed data and their 
uncertainties. The range of values for rjMT in the second col- 
umn of Table 2 is calculated at 2 cr level, i.e. by selecting all 
the simulations for which < Xbost + 4, where Xbest is the 
value for the best fit solution. 

Inspection of Table 2 reveals a clear dichotomy between 
the cases of 47 Tuc, M3 and NGC 6752 and that of lj Cen. 
The 2cr intervals of confidence for t^mt for the first three 

ception was made only for the peculiar case of the velocity dis- 
persion in NGC 6752 (see Section 4.3). 



GCs are in agreement with an almost even distribution of 
MT-BSS and COL-BSS (a slight predominance of COL-BSS 
being indicated by the best fit models). Instead, in the best 
fit model of uj Cen only 14% of BSS are COL-BSS and the 
2 a interval of confidence is compatible with rjuT ~ 1. 

In the following, we will discuss the case of these three 
clusters separately, whereas ui Cen will be examined in Sec- 
tion 5. 



4.1 M3 

In the case of M3, our findings contrast earlier suggestions 
by Sigurdsson et al. (1994). They hypothesized that the BSS 
observed in the peripheral regions of M3 could be COL-BSS 
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born initially in the core and then ejected in the periphery 
because of their natal kick. Indeed, our dynamical simula- 
tions show that COL-BSS are ejected from the entire cluster 
(if the kick is too high) , or sink back to the core in < 1 Gyr, 
if the kick is low enough to retain them in the cluster poten- 
tial well. Hence, they cannot account for the peripheral BSS. 
Instead, MT-BSS have in general much more circular orbits 
than the COL-BSS and can remain in the peripheral regions 
of the GC for a time comparable to the typical lifetime of a 
BSS. 



4.2 47 Tuc 

As for 47 Tuc, we confirm the main outcome from paper I: 
i.e. most of the central BSS of 47 Tuc are born from colli- 
sions and all the peripheral BSS originate from mass-transfer 
in binaries. However, in paper I it was suggested that only 
- 25% of the total BSS are MT-BSS. The upward correc- 
tion of this estimate is due to the fact that we have refined 
our parameter grid with respect to paper I. In paper I we 
considered only cases with a relatively low fraction of MT- 
BSS ( ^ 30%) and one single case with rjuT = 1. For this 
paper we explored a more complete parameter space, the 
percentage of MT-BSS (with respect to the total number of 
BSS) going from to 100%. Our present best fit (46% MT- 
BSS and 54% COL-BSS) is still consistent (at 2 a) with the 
resuh of paper I (25% MT-BSS and 75% COL-BSS), given 
the observational uncertainties, especially in the central bins 
(see Fig. 2). 

4.3 NGC 6752 

The parameters of the best fit model for NGC 6752 are 
similar to those inferred for M3 and 47 Tuc. However, for 
this cluster there may be a caveat about the central veloc- 
ity dispersion, whose value is still very uncertain. From the 
integrated-light spectra of the core, Dubath et al. (1997) 
derived a = 4.91^ 4 km s~^. More recently, from proper 
motion measurements of a sample of ~ 1000 stars with 
the Wide Field Planetary Camera 2/ Hubble Space Tele- 
scope (WFPC2/HST), Drukier et al. (2003) obtained a = 
12.4 ± 0.5 km s'V Thus, we ran simulations for each of 
these two possibilities. 

If we assume tiast ~ 1.5 Gyr, the best fit model ob- 
tained for a = 4.9 km s~^ (Fig. 3) is not satisfactory, since 
we cannot reproduce the increasing number of BSS located 
in the cluster outskirts. In fact, a low velocity dispersion 
implies that the dynamical friction time-scale is too short 
(<^ 1 Gyr) and PBs rapidly shrink towards the core. If we 
adopt a = 12.4 km s~^, the rising trend of the BSS radial 
distribution in the cluster periphery can be reproduced; but 
the (~ 5) is still unacceptable, because of the disagree- 
ment with the data points in the inner regions. Since none of 
the two options was satisfying, we ran new simulations for 
other values of a (ranging between 4.9 and 12.4 km s~^), 
until we found an acceptable fit to the observed BSS distri- 
bution. That occurs for a = 7lf km s~^ (x^ = 1-7, reported 
in Table 2), marginally consistent with the measurement of 
Dubath et al. (1997). 

However, for this cluster other marginally acceptable so- 
lutions for the BSS radial profile exist: e.g. assuming fiast = 4 




10 20 30 

Figure 3. Distribution of BSS normalized to the distribution of 
HB stars in NGC 6752. The solid line indicates the observations 
(filled circles; Sabbi et al. 2004). The dotted line (connecting open 
squares) reports the best fit numerical simulation for a = 4.9 km 
s^^; the long-dashed line (connecting open circles) is for a = 7 
km s^^ (reported also in Fig. 2); the dot-dashed line (connect- 
ing crosses) is for the case a = 12.4 km s~^. In all these cases 
tlast = 1-5 Gyr. The short-dashed line connecting open triangles 
is obtained for a = 12.4 km and tiast = 4 Gyr. 

Gyr, a fit {x^ ^ 3) can be obtained also for a = 12.4 km s"'^ 
(open triangles and short-dashed line in Fig. 3). Considering 
this uncertainty and the peculiarities of NGC 6752 (e.g. the 
unusual density profile), the aforementioned best fit value 
a — 71^ km s~^ cannot be taken as an indirect measure- 
ment of the central velocity dispersion in the cluster. A new 
measurement of a is strongly needed to understand the dy- 
namics and evolution of NGC 6752. 

4.4 The location of the minimum in the BSS 
radial distributions 

The radial distribution of BSS (normalized to the number 
of RGB or HB stars) in M3, 47 Tuc, and NGC 6752 clearly 
displays a minimum at a distance from the cluster center in 
the range 5 — 10 rc. The area surrounding this minimum was 
dubbed 'zone of avoidance' (paper I). We expect that its lo- 
cation roughly corresponds to the radius rav below which all 
the PBs of mass > mess (where mass is the minimum mass 
for generating BSS via mass transfer) have already sunk to- 
wards the cluster center because of dynamical friction. In 
other words, there is a small probability of observing a MT- 
BSS located at a distance from the cluster center of the 
order of r-av. Most likely, the MT-BSS which were originally 
located at r < rav can now be detected at r < rc. 

The time of dynamical friction tdf for a mass mess lo- 
cated at a radius r from the center of a GC is given by 
(Binney & Tremaine 1987) 



4 lnAG2(2 7r)i/2 mBSSp(r)' ^' 
where In A ~ 10 is the Coulomb logarithm and G the gravi- 
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Table 2. Mass Transfer vs CoUisional Scenario 





'?MT 


Vmt{< rc) 


r?MT(> l-av) 




''^av /'^c 


NGC 104 
NGC 5272 
NGC 6752 
NGC 5139 


46+°-°^ 

Q 4^^+0.16 
"■^-^-0.20 


r| + 0.002 
r|+0.002 

02+" °^ 
'^■'-'^-0.04 


Q 91+0.02 
07+0.01 

Q9Q+0.10 


0.5 
2.2 
1.7 
0.8 


9.0 
4.9 

12.9 

1.05 



The second column reports the fraction rjMT of MT-BSS with 
respect to the total number of BSS (the errors are at 2 cr, i.e. 
referred to = '^bost '^)' '^^^ third (fourth) column reports 
the fraction of MT-BSS inside rc (outside rav) with respect to the 
total number of BSS inside rc (outside rav). The fifth column lists 
the reduced obtained for the best fit models. The rightmost 
column reports rav normalized to rc. 
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tational constant. p(r) and a{r) are the cluster density and 
the velocity dispersion at the radius r. We note the depen- 
dence of tdi on the density p(r) and especially on the cube 
of the velocity dispersion. Given the definition of this time- 
scale, Tav can be calculated setting t^f = fgc, where t^c is the 
hfetime of the cluster (~ 12 Gyr) , and using the simulated 
stellar density profile and the velocity dispersion distribu- 
tion of the cluster for inferring rav from p(rav) and a (rav). 
In doing the calculation, we choose mess = 1.2 M©. 

The results (for rav normalized to rc) are shown in the 
last column of Table 2 and in Fig. 2 (vertical dotted line). 
The values of rav are always consistent with the zone of 
avoidance observed in the data. This confirms for M3 and 
NGC 6752 what had been found for 47 Tuc (paper I). 

Third and fourth columns of Table 2 report the values 
of the fraction rj-^^rj, (i) for the MT-BSS located within the 
GC core rc (r7MT(< r^)), and (ii) for the MT-BSS found at 
a distance from the cluster center larger than rav (?7mt(> 
J'av)). These findings indicate that the position of a BSS in 
a GC (at least in those showing a bimodal distribution of 
the BSS) can be used as a strong indication of the nature of 
the progenitor. In particular, a BSS which is located outside 
the zone of avoidance almost certainly is a MT-MSS; on the 
contrary, a BSS found close to the cluster core is a COL-BSS. 

Among the considered clusters, only u Cen does not 
show any zone of avoidance. If we calculate the expected 
radius of avoidance also for this cluster, we find rav ~ r^. 
This indicates that in uj Cen the dynamical friction is far 
less efficient than in the other clusters. In uj Cen the value of 
^mt{> '"av) is similar to the other clusters, whereas 771^^^. 

(< 

rc) is much higher (see Section 5). 

We note that very similar values of both t;mt and rav/''c 
have been obtained for an intermediate density GC (M3) 
and for two much denser clusters (47 Tuc and NGC 6752). 
The analogy between these three clusters is even more ev- 
ident in Fig. 4, where radial distances are normalized to 
the avoidance radius, rav The radial distributions of BSS in 
M3, 47 Tuc and NGC 6752 nearly superimpose. This also 
indicates that rav is a crucial parameter in describing the 
bimodal distribution of BSS. 



Figure 4. Observed distribution of BSS normalized to the dis- 
tribution of HB (47 Tuc, NGC 6752) or RGB stars (M3). On the 
X-axis we plot the distance from the center normalized to rav 
Crosses indicate the observational data for 47 Tuc, open squares 
for M3 and open triangles for NGC 6752. The observational points 
are connected by a dashed line for M3, a dotted line for 47 Tuc 
and a solid line for NGC 6752. 



5 THE CASE OF uo CENTAURI 

The main characteristic of the BSS radial distribution of 
Lu Cen, with respect to the other clusters discussed here, is 
the absence of any central peak. In Fig. 2 we show the high- 
resolution portion of the data set presented by Ferraro et 
al. (2006a): the distribution appears to be flat. This is quite 
unusual, since in all the other GCs surveyed up to now the 
BSS distribution is peaked at the cluster center. Fig. 6 by 
Ferraro et al. (2006a) suggests that this behaviour extends 
up to at least ~ 7rc (~ 20' from the cluster center). 

This peculiarity adds to the many other unique features 
displayed by lo Cen (see Section 3); but our simple model 
provides a satisfactory fit to the data (see Fig. 2 and Table 2) 
also in this case. However, at variance with the results for 
the three clusters examined in Section 4, the best fit model 
requires that a large majority of the BSS are born from PBs 
(see Table 2). In other words, only a mere ~ 14% of the 
BSS detected in ui Cen may have a coUisional origin, i.e. we 
expect that only ~ 44 of the 313 BSS reported by Ferraro 
et al. (2006a) are COL-BSS. 

The alternative hypothesis that all the others BSS 
formed via collisions have been ejected from u Cen sounds 
unrealistic, since the central escape velocity of u) Cen (~ 40 
km s~^) is comparable with that of other clusters in our 
sample and there is no motivation for assuming that colli- 
sions in LU Cen should generate larger recoil velocities with 
respect to other clusters. As a consequence, the dynamical 
modelling of the BSS presented in this work predicts that 
COL-BSS are now produced with a low rate in ui Cen. 

Are there viable explanations for this underproduction 
of COL-BSS in u Cen? We can just note that the charac- 
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teristic radius rav for ui Cen is much smaller (both in terms 
of rc and of physical units) than in the other clusters of our 
sample (rav ~ instead of the typical value ranging be- 
tween 5 Tc and 15 rc). This means that only a small number 
of PBs had enough time to sink into the core due to mass 
segregation. 

PBs that sink within the core are thought to be the 
progenitors (as far as modified by three body interactions) 
of most of the binaries hosted in the core of current GCs, 

since the formation of binaries from two-body interactions is 
a quite unlikely process. Because of their large cross-section, 
core binaries suffer repeated three- and four-body interac- 
tions, eventually exchanging companions, being ionized or 
hardened. Thus, core binaries are required to form COL- 
BSS, via three- or four-body encounters. Then, the relative 
lack of binaries in the core of ui Cen, due to the inefficiency 
of dynamical friction, could have quenched the formation of 
COL-BSS. 

This view could be strongly supported by the observa- 
tion of a low (<C 10 %) fraction of binaries in the core of 
u) Cen. Unfortunately, this fraction cannot be derived with 
classical methods (Rubenstein & Bailyn 1997; Bellazzini et 
al. 2002), because of the wide spread in rnetallicity of the 
stars in this GC (Norris, Freeman & Mighell 1996; Suntzeff 
& Kraft 1996). An indirect indication of the occurrence of 
a low fraction of binaries in the core of u} Cen comes from 
inspection of Fig. 2, which shows the radial distribution of 
the X-ray sources observed in cu Cen by XMM (Gendre, Bar- 
ret & Webb 2003), normalized to the number of RGB stars. 
Most of the X-ray sources are believed to be cataclysmic 
variables or low mass X-ray binaries (Gendre et al. 2003). 
So they should be among the most massive objects in the 
GC. Their flat radial distribution further supports the hy- 
pothesis that mass segregation in lo Cen has not driven yet 
a sizeable number of massive binaries in the central region. 



6 NUMBER OF BSS VERSUS My 

Until now, we only considered the relative frequency of BSS, 
i.e. the number of BSS normalized to that of HB and/or 
RGB stars. But what happens if we consider just the number 
of BSS hosted in each cluster, without any normalization? 

First of all, it is worth noting that NGC 6752 {lo Cen) 
hosts about three times fewer (more) BSS than either 47 Tuc 
or M3 (Fig. 5). We can ask whether this is at all correlated 
to the mass of the host cluster, or, equivalently, to its abso- 
lute magnitude My'? Plotting in Fig. 5 the total number of 
observed BSS per cluster against My for our four clusters, 
we do not see any straight correlation. This result agrees 
with the fact that Piotto et al. (2004) found evidence, in a 
large sample of GCs, of a lack of correlation of the number of 
observed BSS with either stellar collision rate (as would be 
expected if all BSS were collisional in origin) or mass (which 
is proportional to the collision rate). 

This fact led Davies et al. (2004) to propose that BSS 
are made through both channels, suggesting that the number 
of COL-BSS tends to increase with cluster mass, while MT- 
BSS tend to decrease with total mass, as PBs drifting by 
dynamical friction in the core had already time to burn as 
BSS in the past. The combination of these two competing 
behaviours was found to produce a population of BSS weakly 



dependent on the total mass and core collision rate of their 
host GC. 

Fig. 5 sketches the predictions by Davies et al. (2004). 
In particular, the dashed (dotted) line indicates the contri- 
bution expected from the MT- (COL-)BSS; the solid line is 
the total. 

47 Tuc and M3 lie close to the solid line and, according 
to the Davies et al.'s model, their current BSS population is 
obtained by a blending of COL-BSS and MT-BSS contribut- 
ing in roughly equal rmrnber. This fully agrees with what we 
found firom the comparison between our simulations and the 
radial distribution of BSS. 

In contrast, our model and Davies et al.'s scenario dis- 
agree in the case of both NGC 6752 and lu Cen. NGC 6752 
lies well below the solid curve. According to Davies et al. 
(2004), its BSS should predominantly be MT-BSS; whereas 
in our model NGC 6752 should host MT-BSS and COL-BSS 
nearly in equal number. In the Davies et al.'s model a; Cen is 
expected to have only COL-BSS, while we have shown that 
MT-BSS should dominated 

Moreover, though the picture proposed by Davies et al. 
(2004) is interesting, detailed cluster-to-cluster comparisons 
suggest a much more complex scenario, where the dynami- 
cal history, the original PB content and the current dynam- 
ical state of each cluster seem to play a major role (Fer- 
raro et al. 2003a). Indeed, clusters with the same integrated 
magnitude harbor quite different BSS populations. In par- 
ticular, two GC pairs offer the possibility of demonstrating 
the complexity of the emerging scenario: (i) M3 and M13 
are almost twins GCs (Ferraxo et al. 1997b). They have the 
same integrated magnitude {My ~ —8.6), same metallicity 
([Fe/H] — —1.6), same mass, but display a quite different 
BSS content (Ferraro et al. 2003a). In particular, M13 har- 
bors a factor 5 fewer BSS than M3. (ii) NGC 6752 has an 
integrated magnitude (My ~ —7.7) which is quite similar 
to M80 (My ~ -7.9); but again the BSS content in the two 
clusters turns out to be quite different: NGC 6752 harbors a 
BSS population which is a factor 10 lower than that found in 
M80 (note that the BSS population of M80 is comparable in 
size to that found in uj Cen). Ferraro et al. (1999) suggested 
that the anomalously large population of BSS in M80 could 
have originated in the core-collapsing phase of this cluster. 

The model discussed in this paper is based on a differ- 
ent approach from that of Davies et al. (2004). Neither this 
work nor Davies et al.'s accounts for the dynamical evolu- 
tion of the cluster; but we consider additional information 
in comparison with Davies et al. (2004), i.e. the observed 
BSS radial distribution. 

Davies et al. (2004) try to predict the properties of the 
BSS from the present properties of the host cluster, and 
these properties might not reflect all the stages in the clus- 
ter evolution. In particular, they risk overlooking the im- 
portance of the earlier evolutionary phases which can have 
a strong impact on the characteristics of BSS population. 

In our approach the cluster evolution is frozen for the 
last 2 Gyr only, when the cluster properties are not expected 
to have changed significantly. All the effects of the previous 
evolution of the cluster onto BSS are intrinsically stored in 



^ Davies et al. (2004) admit that their model is not adaptable to 
slow evolving GCs like ui Cen. 
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Figure 5. Number of BSS as a function of the cluster ab- 
solute magnitude My. Open triangle refers to NGC 6752, 
filled triangle to M80, open square to M3, filled square to 
M13, cross to 47 Tuc, filled circle to lo Cen. The val- 
ues of My come from the Harris catalogue (Harris 1996; 
http : //physwww.mcmaster . ca/'/,7Eh£irris/mwgc.dat). The solid 
line shows the model by Davies et al. (2004), while the dashed 
(dotted) line refers to the MT-BSS (COL-BSS) according to the 
same model. 

the initial parameters which we impose on the BSS popula- 
tion: masses, lifetimes, velocities, locations and the amount 
of injected MT-BSS and COL-BSS. Hence, we can use the 
observed BSS radial profile for inferring the ratio between 
the injected MT-BSS and COL-BSS (which is related to 
their location), and for investigating the effects of the clus- 
ter (present) properties on that (e.g. what is the importance 
of collisions in the recent life of the cluster; and, what is the 
residual fraction of PBs needed to remain in periphery up 
to now). 

Both these methods need refinements; however, the ap- 
proach presented in this paper has the advantage that it is 
able to fit the BSS population also in clusters where Davies 
et al.'s method appears inaccurate. A complete understand- 
ing of the relation between a cluster and its BSS population 
will likely be possible only by running N-body simulations 
of the cluster from its formation, accounting also for three- 
body encounters, PB evolution, stellar evolution, etc. 



the central region of the cluster. The MT-BSS are slightly 
less abundant than the COL-BSS, but populate all the GC. 
The density of BSS reaches a minimum in a so-called zone 
of avoidance, which separates the portion of the GC mostly 
occupied by COL-BSS from the cluster outskirts, where the 
MT-BSS dominate. The location of the zone of avoidance is 
explained by accounting for the effects of the dynamical fric- 
tion on the PBs which were massive enough for generating 
the observed BSS. 

The picture described above can also be applied to 
Cen; but in this case the lack of a central peak in the 
BSS radial distribution requires that the large majority of 
the BSS derive from PBs. The very low rate of production of 
COL-BSS could be in turn attributed to the fact that mass 
segregation has not yet driven a sizeable number of PBs to 
the central region of the cluster to produce BSS. 

A very interesting further development of this research 
will be to perform a comparison between the location of a 
significant sample of BSS in a GC and their spectroscopic 
properties. According to the findings of this work, the posi- 
tion in the GC might represent a strong dynamical clue for 
the formation mechanism of a given BSS. If it is located out- 
side the zone of avoidance, the BSS almost certainly results 
from evolution of a PB; if it is harbored in the cluster core, 
the BSS has most likely a coUisional origin. On the other 
hand, indication about the origin of the same BSS can be 
independently obtained from high resolution spectroscopy. 
Indeed the chemical signature of the MT-BSS formation pro- 
cess has been recently discovered in 47 Tuc (Ferraro et al. 
2006b). The acquisition of similar sets of data in clusters 
with different structural parameters and/or in different re- 
gions of the same cluster will provide an unprecedented tool 
for confirming the scenario presented here and to finally ad- 
dress the BSS formation processes and their complex inter- 
play with the dynamical evolution of the cluster. 



8 ACKNOWLEDGMENTS 

We thank Craig Heinke, Fred Rasio and Emanuele Ripa- 
monti for useful discussions. We acknowledge the anonymous 
referee for the critical reading of the manuscript. M. M. ac- 
knowledges the Northwestern University and the Center for 
Gravitational Wave Physics (Pennsylvania State University) 
for kind hospitality. M. C. and A. P. acknowledge grant 
PR1N05 2005024090-002. This research was supported by 
the contract ASLINAF 1/023/05/0. 



7 SUMMARY 

We have exploited the recent determination of the radial dis- 
tribution of BSS in four GCs, in order to investigate which 
mechanism of BSS formation prevails in these stellar sys- 
tems. Our conclusion is that the two main formation paths 
proposed so far, i.e. mass-transfer in PBs and merging of MS 
stars due to collisions in the cluster core, must co-exist and 
have similar efficiency both in a low density cluster (MS) 
and in much denser clusters, like J^l Tuc and NGC 6752. 

In particular, in M3, 47 Tuc, and NGC 6752 the COL- 
BSS sum to ~ 50 — 60% of the total and mostly reside in 
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